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» ABSTRACT 


Expressions  for  a generalized  Modified 
Transitional  Butt 9r worth-Chebyshev  (MTBC)  filter  are 
derived.  The  characteristics  of  this  filter  as  applied 
to  digital  filter  design  are  investigated.  It  is  shown 
that  by  adjusting  location  and  order  of  the  inserted 
zeros,  the  cut-off  slope  rate  of  the  filter  can  be 
traded  for  maximum  attenuation  in  the  stop-band.  | 

The  performance  of  this  MTBC  filter  is  compared  to  I 

i 

that  of  Butterworth,  Chebyshev,  transitional  j 

Butterworth-Chebyshev  filters  together  with  those  ' 

suggested  by  other  investigators  [1]-[3].  It  is  shown 
that  the  stop-band  attenuation  can  be  significantly 

I 

increased  without  great  sacrifice  of  cut-off  slope  i 

I 

rate. 

Step  response  of  this  MTBC  filter  is  also  obtained 
and  compared  with  other  filters.  Various  tabulations  | 

as  well  as  graphs  of  this  filter  are  given  for  design  t 

purposes.  A computer  program  is  developed  for  the 
design  of  this  filter. 
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INTSODOCTION 


The  digital  filter  is,  as  defined  by  Rabiner  et  al  [10], 
'•  a computational  process  or  algorithm  by  which  a digital 
signal  or  sequence  of  numbers  (acting  as  input)  is 
transformed  into  a second  sequence  of  numbers  termed  the 
output  digital  signal”. 

The  area  of  digital  filtering  can  be  divided  into  two 
major  subdivisions  as  Finite  Impulse  Response  (FIR)  filters 
and  Infinite  Impulse  Response  (HR)  filters. 

During  the  development  of  digital  signal  processing,  the 
interest  cf  the  investigators  in  HR  and  FIE  filters  varied. 
Before  the  introduction  of  the  FFT  algorithm  by  Cooley  and 
Tukey  ( 1965)  HR  filters  were  much  more  efficient  than  FIR 
filters.  Stccham's  work  [13]  on  the  FFT  method  of 
performing  convolution  indicated  that  implementaticn  of 
high-order  FIR  filters  could  be  made  extremely 
computationally  efficient  ; thus,  comparison  between  FIR  and 
HE  filters  are  no  longer  strongly  biased  toward  the  latter 
[5].  Beoause  FIR  filters  require  very  high  orders  to 
produce  a sharp  attenuation  shape,  they  are  not  often  used 
for  real-time  filtering  of  waveforms.  Recently,  due  to  the 
increase  in  computing  capabilities  in  digital  signal 
processing  and  the  availibity  of  long  charge  transfer  device 
(CTD)  tapped  delay  lines  (TDL) , FIR  filters  are  favored  over 
HE  filters.  However,  in  applications  like  design  of 
digital  comb  filters  HR  filters  are  the  unique  alternative. 

There  are  three  basic  design  techniques  of  HR  digital 
filters  [5]. 
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which  is, 
approximate 


First  method  is  the  direct  design, 
appropriately  placing  poles  and  zeros  to 
required  frequency  response. 

A second  method  is  to  use  an  optimization  procedure  to 
place  the  poles  and  zeros  to  match  arbitrary  frequency 
response  specifications. 

Finally  the  third  technique  makes  use  of  highly  advanced 
art  of  continuous  filter  design.  This  technique  of 
designing  digital  filters  from  continuous  filters  by  means 
of  mathematical  transformations  is  the  most  popular  HR 
digital  filter  design  technique. 

Standard  Z-transform, Bilinear  Z-transform  and  the 
matched  Z-transform  make  possible  direct  transformation  from 
S-domain  to  2-domain,  preserving  essential  characteristics 
of  analog  frequency  response. 

Existence  of  frequency  transformations  reduces  tae 
problem  to  design  a frequency  normalized  prototype  low-pass 
filter. Then  using  appropriate  frequency  transformation,  this 
prototype  nay  be  converted  into  desired  band-pass, 
band-reject  or  high-pass  filter.  Popular  prototype  filters 
are  Butterworth,  Chebyshev,  elliptic  and  hybrid  transitional 
filters.  Frequency  transformations  for  digital  filters  are 
discussed  in  various  literatures  ( [6]  and  [7]  ). 

The  problem  of  designing  low-pass  prototype  filters, 
which  possesses  better  stop-band  attenuation  and  cut-off 
slope  characteristics  than  existing  prototypes  has  always 
attracted  the  researchers  in  the  signal  processing  area. 

Budak  and  Aronhime  suggested  [1]  modification  of 
maximally  flat  rational  functions  by  introducing  a pair  of 
finite  transmission  zeros  such  that  the  maximally  flat 
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characteristic  is  maintained  but  tne  cut-off  slope  can  be 
made  steeper  without  great  sacrifice  of  stop-band 
attenuation. 

Dutto  Roy  [2]  investigated  a more  general  case  allowing 
insertion  of  multiple  pairs  of  transmission  zeros,  either 
coincident  or  distinct. 

Introducing  multiple  pairs  of  jw-axis  zeros  in  all  pole 
Chebyshev  transfer  functions  are  investigated  by  Agarwal  and 
Sedra  [ 3 ]. 

The  most  attractive  feature  of  these  finite  zero  filters 
is  that  they  offer  the  filter  designer  a great  degree  of 
freedom  in  choosing  the  location  and  order  of  the  zeros  to 
trade  cut-off  slope  for  stop-band  attenuation. 

In  thiiS  thesis,  a modified  Transitional 
Butterworth-Chebyshev  filter  is  developed,  which  is  a more 
general  case,  introducing  finite  coincident  or  distinct 
multiple  pairs  of  transmission  zeros  in  transitional 
Butterworth-Chebyshev  filter. 

Trade-off's  between  tne  order  of  the  filter,  the  order 
of  transmission  zeros,  stop-band  attenuation  and  cut-off 
slope  are  pointed  out.  Graphs  helpful  in  the  design  of  such 
filters  are  obtained. 

Performances  of  Butterworth,  Chebyshev,  Transitional 
Butterworth-Chebyshev  filters  and  the  designs  suggested  in 
references  [2]  and  [3]  are  compared  with  those  represented 
in  this  thesis  for  the  orders  of  three  through  eleven.  A 
computer  program  is  developed  to  implement  the  filters 
mentioned  above. 

In  addition,  the  time-domain  response  of  digital  filters 
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is  studiad.  There  are  many  applications,  such  as  digital 
HTI  filters,  for  which  one  is  interested  in  the  transient 
responses  cf  filters  that  are  specified  in  the  frequency 
domain.  Step  responses  of  the  filters  that  are  discussed  in 
this  report  are  plotted  and  compared. 
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II 


• DERIVATION  OF  MODIFIED  TRANSITIONAL 
BDTTERWORTH-CHSBYSHEV  FILTERS 


A.  INTRODUCTION 


The  most  popular  technique  for  designing  HR  digital 
filters  is  to  digitize  an  analog  filter  that  sarisfies  the 
design  specifications  [5],  There  are  many  techniques  for 
designing  analog  low-pass  prototype  filters.  Among  the  well 
known  analog  filter  classes  are  the  maximally  flat 
(Eutterworth)  and  equal  ripple  (Chebyshev)  filters. 

Butterworth  filters  are  simple,  excellent  in  the 
pass-band  and  monotonic  in  both  pass-band  and  stop-band.  The 
Chebyshev  filters  are  superior  at  and  near  cut-off  frequency 
and  at  stop-band. 

The  transitional  Butterworth-Chebyshev  (TBC)  filters 
combine  the  desirable  attributes  of  these  two  filters  in  a 
single  approximation  that  is  given  by 


I 


where 


k=  Weighting  factor 
n=  Order  of  filter 

C ,(w)=  (n-k)"^  order  chebyshev  polynomial. 

When  k=n,  |F(jw)  |2  is  identical  with  the  Butterworth 
function.  When  k=0,  |F(jw)|2  is  identical  with  the  Chebyshev 


function.  With  a varying  value  of  k,  a TBC  filter  possesses 
sene  characteristics  of  each.  As  k approaches  to  n,  a TBC 
filter  behaves  more  like  a Butterworth  filter,  as  k 
approaches  to  zero,  it  behaves  more  like  a Chebyshev  filter. 

In  this  chapter,  modification  of  TBC  filters  by 
introducing  finite  coincident  or  distinct  multiple  pairs  of 
transmission  zeros  will  be  discussed.  It  will  be  shown 
that,  using  a weighting  factor  k,  and  the  location  and  order 
of  inserted  zeros  as  parameters,  attenuation  in  the 
stop-band  may  be  traded  for  sharpness  of  the  cut-off 
characteristics. 

Expressions  for  the  cut-off  slope  and  minimum 
attenuation  in  the  stop-band  are  derived  in  terms  of  order 
of  the  filter,  weighting  factor,  location,  and  order  of 
inserted  zeros. 


B.  MODIFICATION  OF  THE  TBC  FILTEHS  WITH  COINCIDENT 
TRANSMISSION  ZEROS 


Introducing  m identical  pairs  of  transmission  zeros  at 
±jWo  to  eg.  (1-1),  we  have 


In  order  to  normalize  |F(jw)  |^,  i.e,  to  force  |F(jw)  to  be 
egual  to  1/2  at  w= 1 , the  constant  K should  be 
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Then  eq.  (2-1)  becomes 


we 


Taking  the  derivative  of  eg.  (2-10)  with  respect 


to  w. 


At  the  cut-off  frequency,  w=1,  we  have 


f '(lo)  ■= (l-ii) 

K(u))=  i 

k7ui')=  (l-w) 

|F(j>o}|  = M 


Substituting 
we  obtain 


these  values  into  eg.  (2-11)  and  simpiif lying. 


rvi  k.4  (n-  It')  ^ 

^ (u)o^-0 


For  k=n , the  result  agrees  with  the  cut-off  slope  of  MB 
function,  which  is  derived  in  [2].  For  k=0,  the  result 
agrees  with  the  cut-off  slope  of  MC  function  as  derived  in 
[3]. 


3 • Stop^band  characteristics 

In  general,  stop-band  characteristics  of  finite  zero 
filters  will  te  of  the  fora  shown  in  FIG.1^. 


♦ Because  of  the  inherent  limitations  in  the  plotting 
subroutine  utilized  to  provide  the  graphs  in  this  thesis,  it 
was  necessary  to  add  supplementary  axes  to  show  proper 
scaling  for  seme  of  the  graphs. 
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The  zeros  introduced  will  cause  a peak  in  the  stop-band,  at 
a frequency  w^>w, . The  minumum  attenuation  in  the  stop-band 
may  he  defined  as 


MT8C  ™ T p 


where 

In  the  stop-hand  eg.  (2-4)  becomes 


) = 


f (ijf)  = (iOo'-U)>) 

w(io,)= 


(2.2I) 


Taking  derivatives  of  these  values  and  substituting  them 
into  eg.  (2-11)  results  in 

Ip,.  ■ ] 

(2-23) 


Combining  eg.  (2-19)  and  eg.  (2-23),  we  obtain 


n-  imrt 


(2-2^ 


Thus 


(too'-O""'  ^^(-'-"'-'■-0 


i(n-2rt\') 


,f\-2m  2m 
*T1 


I 


■ 


F 


And  minmnuin  stop-band  attenuation  will  be  given  by 

C^MTac  **  10  n Ic^  C n 2 ^ (^o) 
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Plots  of  stop-band  attenuation  and  cut-off  slope  of 
MT2C  filters  with  two  coincident  transmission  zeros,  for 
orders  3 thrcugh  11  are  given  in  figures  2 and  3. 


C.  MODIPICAIION  OF  T3C  FILT2SS  WITH  DISTINCT  THANSMISSION 
ZEEOS 


-+W 

Consider  the  n—  order  TBC  function  with  m pairs  of 
finite  distinct  zeros  at±jw;  , where  i=1,2, ,a  and  w;  >1. 

The  iragnitude  sguared  function  of  TBC  function  w^th 
distinct  tr ansmission  zeros  will  be  given  by 


m 


r(j^) 


.71 


(“!) 


Let 


i»1 


3tu>)=  J 

< £ I 

Putting  these  values  into  eg.  {2.-21)  results  in 
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f(jw) 


9Cu)) 

k(u>) 


(..i) 
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'•  • stop-band  and  cu^of | chacacteriot ics 


occur,  is  given  by 


iFCju.)!' 


-=  0 


(d  a-vdii 


respect  tc  w,  we  obtain 

iZCn-k-V^  2.n-\ 

I u f I . \ 1 f/ 

F(ju.) 


_ 3(10) 


KVio)[i-t  iC<^) 


solved  for  Wp  may  be  found  to  be 


^ U), 


m 


I 


i- 


+ ri  = 0 


i-band  peaic 

will 

(ii?) 

eg.  (2-28) 

with 

1 V (10)  J 
I'L 

K(u))  i 

H 

equation 

to  be 

(..39 

identical. 

attenuation 

Using 

peaks 

may  be  fcund  to  be 

* Jo  JL03 

where 

= |F(ju)pi) 

and  miniaium  stop-band  atrtenuation  is  given  by 


^ A\i*i  = »v\a7l  1-'] 


H 
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Using  eg.  (2~5)  and  differentiating  eg.  (2-27)  with  respect 
to  w,  the  cut-off  slope  becomes 

which  agrees  with  eg.  (2-17)  when  the  w-,  's  are  identical. 

Plots  of  cut-off  slope  and  stop-band  attenuation  of 
MTBC  filters  with  two  distinct  transmission  zeros,  for 
orders  3 through  11  are  given  in  figures  3-17. 

D.  SUMMAKY 

The  Transitional  B utterworth-Chebyshe v filter  combine 
the  best  features  of  the  Butterworth  and  Chebyshev  filters. 
A Modified  Transitional  Butterworth-Chebyshe v filter, 
obtained  by  introducing  finite  transmission  zeros  to  a 
Transitional  Eutterworth-Chebyshev  filter,  possesses  cut-off 
slope  and  stop-band  attenuation,  which  are  dependent  on  the 
modification  parameters  Wo  and  m.  The  closer  the  Wa  is  to 
unity,  the  steeper  the  cut-off  slope.  However,  this 
improvement  in  the  cut-off  region  results  in  degradation  of 
stop-band  attenuation.  Thus,  using  Wo  and  m as  parameters, 
an  advantageous  trade  between  attenuation  in  stop-band  and 
sharpness  of  the  cut-off  characteristics  can  be  made. 
Graphs  are  given  to  serve  as  guides  in  trading  cut-off  slope 
for  stop-band  attenuation. 
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Figure  5 - COT-OFF  SLOPE  OF  MTBC  FILTER  WITH  TWO  DISTINCT 

ZEROS  (w  =1«06,  n=6) 
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Figure  6 - CUT-OFF  SLOPE  OF  ."ITBC  FILTER  WITH  TWO  DISTINCT 

ZEROS  < w =1*06,  n=7  > 
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Figure  7 - CUT-OFF  SLOPE  OF  MTBC  FILTER  WITH  TWO  DISTINCT 

ZEROS  <w  =1*06>,  n=8> 
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Figure  10  - CUT-OFF  SLOPE  OF  MTBC  FILTER  WITH  OISIINCT 
TRANSMISSION  ZEROS  =1*06  ,n=11> 
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Figure  11  - STOP-flAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <w  =1«06,  a=5> 
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STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 


DISTINCT  TRANSMISSION  ZEROS  (w  =1«06.  n=6 > 


Figure  13  - STOP-BAND  AITBNUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  =1*06,  n=7J 
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Figure  13  - STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  (w  =1*06,  n=7> 


36 


■oti-  C$1-  'Ccl-  "t 


Figure  14  - STOP-BAND  ATTENUATION  OF  ilTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <w  =1*06,  n=8> 
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Figure  15  - STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  TRANSMISSION  ZEROS  <w  =1«06,  n=9 > 


38 


Figure  17  - STOP-BAND  ATTENUATION  OF  MTBC  FILTERS  WITH  TWO 
DISTINCT  THANSNIS3ION  ZEROS  =1*06,  n=  1 1 > 
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III.  COHFARISON  OF  MIBC  FUNCTION  WITH  AND  TBC 

FUNCTIONS 


A.  INTRODUCTION 

In  chapter  2,  derivation  of  MT3C  filter_is  given.  It  is 
shown  that,  stop-band  attenuation  and  cut-off  slope  rate  of 
this  filter  depend  on  modification  parameters  m,  k,  and  Wo . 
So  the  filter  designer  have  the  flexibility  to  trade  cut-off 
slope  for  stop-band  attenuation  by  changing  these 
parameters,  without  changing  the  order  of  the  filter. 

In  this  chapter,  performance  of  MTBC  filter  is  compared 
with  a,  MB,  C,  MC,  and  TBC  filters. 

Formulas  for  stop-band  attenuations  and  cut-off  slopes 
of  all  these  filters  are  given  in  Table  I. 
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TABLE  !•  FORMULAS  FOR  CUT-OFF  SLOPES  AMD  STOP-BAND 

ATTENUATIONS 


TABLE  !•  CONTINUED 


3 


From  eg.  (3-1)  ratio  of  cut-off  slopes  of  these  two 
functions  may  be  found  as 


'MTftC/6 


XuiTftC  rt  , 


v.s.  Wg  for  n=3  through  11  and  ra=1,2  are  given  in 
Figures  18  and  19.  Stop-band  attenuation  of  MTBC  may  be 
written  as 

^MT8c  = c.^8+  6(n-rvi-le-i)  4 o20  ity\  Log  ~ (s-3^ 

Then  the  difference  between  stop-band  attenuations  is  given 
by 

■ 

^M-rftc/a=  ^ 6 ■=  ^ (^n- m- k-l) fvi  JLllIL. 

— 


Plots  of  d. 


t/a  Wo  for  n=3  through  11  and  m=1,2 


are  given  in  FIG  20  and  21. 

Figures  19-21  indicate  that,  MTBC  filter  can  be  made  10 
times  steeper  than  B filter  still  having  40  dB  more 
attenuaticp  at  the  stop- band. 


C.  MTBC  V.S.  MB  FUNCTION 


Using  Table  I,  cut-  off  slope  ratios  of  MTBC  and  MB 
filters  may  be  found  as 


_ J2m+  (u3o'*--0  r_Crt-fc)  S k] 


Piets  of  r, 


flrVl  +rv  (^u)o’’-4) 

’-JmA  *'o  for  n»3  through  11  and  m=1,2 


given  in  figures  22  and  23. 

The  difference  between  the  stop-band  attenuations  of  MTBC 


t 


r 


and  MB  is  given  by 

<^MT6c/uft  ^ 

Equation  (3-7)  indicates  that,  stop-band  attenuation 
difference  between  these  filters  doesn't  depend  on  Wo . For 
n*10,  a=  1 , k=1,  with  48  dB  aore  attenuation,  cut-off  slope 
ratio  may  be  changed  from  2.5  to  8.0  by  changing  Wo  from 
1.06  to  2.4. 


D.  MTBC  FUNCTION  V.S.  C FUNCTION 


The  ratio  of  cut-off  slope  of  MTBC  function  to  cut-off 
slope  of  C function  is  given  by 


'HuTic/c 


i + 


jZm 

(ido^'O 


It,  0e+l-.2rt) 


Plot  of  V.S.  Wo  for  n=3  through  11  and  n=1,2  are 

given  in  Figures  24  and  25.  The  difference  between 
stop-band  attenuations  of  MTBC  and  C functions  is  given  by 


d 


MTBC/c 


= -+-,20nh  1-0^ 


ri  Uo^  J 


Plots  of  V.S  Wo  for  n=3  through  11  and  n=1,2  are 
given  in  Figures  26  and  27.  The  Chebyshev  filter  is  known 
to  provide  much  steeper  cut-off  slope  than  the  corresponding 
B and  T3C  filters.  Figure  25  shows  that  the  cut-off  slope 
of  the  MTBC  filter  can  be  made  1.4  times  steeper  than  that 
of  the  Chebyshev  filter,  with  n=5,  m=2,  and  Wo=1.06. 


E.  MTBC  FUNCTION  V.S.  HC  FUNCTION 
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The  ratio  of  cut-off  slope  of  MTBC  filter  to  cut-off 
slope  of  ac  filter  is  given  by 


Cu)ci^-4) 

iZ.*v\  Jt  a’-  (uo’--  ( 


Plots  of  v.s.  w,  for  n=3  through  11  and  m=1,2  are 

given  in  Figures  23  and  29. 

The  difference  between  stop-band  attenuations  of  alBC  and  MC 
functions  is  given  by 

(2-'0 

For  a given  k,  <^MTac/MC  allways  constant. 


F.  MTBC  V.S.  TBC  FUNCTIONS 


The  ratio  of  the  cut-off  slopes  of  dTBC  and  TBC 
functions  is  given  by 


*^T6c/T6C  * 1 -1- 


.m 


(u).'-0  [nM 


(j-ll) 


Plots  of  ^ktuc/tzc  v.s.  w.  for  n=3  through  11  and 
given  in  Figures  30  and  31. 


1 = 1 ,2  are 


The  difference  between  the  stop-band  attenuations  of 
these  two  functions  is  given  by 


ci«7ac/Tftc  = - 6m  m loj 


I 


(3-lp 


Plots  of  <^m7*c/t8c  v*s.  w,  for  n=3  tnrough  11  and  in=1,2 
are  given  in  Figures  32  and  33. 
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Dsing  the  location  and  the  order  of  the  inserted  zeros 
and  the  weighting  factor  as  parameters,  the  characteristic 
curve  of  the  MTflC  filter  can  be  made  steeper  than  that  of 
the  conventional  all-pole  filters  without  greatly 
sacrificing  either  stop-band  attenuation  or  flatness  in  the 
pass  band.  A Modified  Chebyshev  filter  provides  slightly 
better  performance  than  the  MT3C  filter,  at  the  expense  of 
substantial  degradation  of  the  pass-band  flatness.  Graphs 
are  given  to  help  in  the  comparison  and  in  determining  the 
numerical  advantages  gained  by  increased  complexity. 
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Figure  19  - RATIO  OF  COT-OFF  SLOPES  OF  MTEC  AMD  3 FUNCTIONS 
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Figure  22  - 


3ATI0  OF  CUT-OFF  SLOPES  OF  MTBC  AND  MB 
FUNCTIONS  V«S*  M <M=1> 
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Figure  25  - RATIO  OF  THE  CUT-OFF  SLOPES  OF  HTBC  AND  C 
FUNCTIONS  V*S«  H <M=2) 
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Figure  32  - THE  DIFFERENCE  BETWEEN  THE  STOP-B 
ATTENUATIONS  OF  HTBC  AND  T3C  FUNCTIONS  V»S«  W < 
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IV. 


COMPUTER  PROGRAM 


A.  INTRODUCTION 

in  the  previous  chapters,  expressions  for  MTBC  filter 
are  derived  and  the  perforasance  of  this  filter  is  compared 
ifith  various  filters. 

In  this  chapter,  A computer  program  is  developed  to 
implement  MTEC  filter.  It  is  pointed  out  that  all  five 
filters,  which  are  used  to  compare  with  the  MTBC  filter,  are 
the  special  cases  of  the  MTBC  filter.  Thus  the  program 
developed  in  this  chapter  for  the  implementation  of  the  MTBC 
filter  may  also  be  used  to  implement  any  one  of  these  five 
filters. 

A sample  problem  is  worked  out  to  illustrate  the  use  of 
the  plots  presented  in  earlier  chapters,  the  use  of  the 
computer  program,  and  to  point  out  the  flexibility  offered 
by  the  MTBC  filter  to  the  filter  designer. 

B.  COMPUTER  PROGRAM 


In  the  most  general  form  the  transfer  function  of 
filter  is  given  by  eg.  (2-27)  , which  is  repeated  here 
convenience. 


f iCw’-O 

t*l  i=i 


MTBC 

for 
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where  ; 


r 


w = Location  of  i—  inserted  zero 
m = Order  of  inserted  zeros 
n = Order  of  filter 

k = Weighting  factor  of  Transitional  Butterworth 
Chebyshev  filters 

When  k=n,  and  m=0,  IF(jv)l^  is  identical  to  the 

Butterworth  function.  With  k=n  and  m^O,  |F(jw)  |2  is 
identical  to  the  modified  Butterworth  function  which  is 
discussed  by  BUDAK  and  ROY  [1]  - [2].  When  k.=0  and  m=0, 
|F(jw)|2  is  identical  to  the  Chebyshev  function.  With  k=0 
and  a#0  Modified  Chebyshev  function  may  be  implemented  which 
is  discussed  by  AGAHVAL  and  SED3A  [3]. 

When  n=0  and  k#0,  |F(jw) |z  is  identical  to  the 
Transitional  Butter worth-Chebyshev  filter.  Finally  allowing 
both  k and  m to  vary.  Modified  Transitional 
Butterworth-Chebyshev  filters  may  be  implemented. 

The  program  consists  of  two  main  parts.  In  the  first 
part  the  analog  filter,  which  is  given  by  eq.  (2-27),  is 
implemented  and  the  frequency  response  is  plotted. 

In  the  second  part,  the  analog  transfer  function  F(jw) 
is  first  piedistorted  then  transformed  into  z-domain  by 
algebraic  substitution  method  (using  Bilinear 

z-transf ormation)  to  obtain  H (z)  . Digital  transfer  function 
[ H(z)  J is  then  factored  into  second  order  cascaded  stages. 
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Finally,  frequency  response  curves  ( |F(jw)|  v.s.  w,  20 

I 

log  |F(jw)l  v.s.  Wa  and  F (Jjw)  v.s.  w* ) are  drawn. 

Bilinear  Transformation  is  preferred  over  the  other 
available  algebraic  substitution  methods  (i.e.  Impulse 
invariant.  Matched  z-transf ormation)  in  obtaining  E (z) , 
mainly  for  the  following  reasons  [4]  : 

(1)  it  has  the  property  that  realizable  stable 
continuous  systems  are  mapped  to  realizable  stable  digital 
filters. 


(2)  Wideband  sharp  cut-off  continuous  filters  can  be 
mapped  to  wideband  sharp  cut-off  digital  filters  without  the 
aliasing  in  the  frequency  response. 


(3)  After  Bilinear  Transformation,  the  relation  between 
the  analog  and  digital  frequencies  is  given  by 


By  choosing  cut-off  frequencies  approaching  to  fj/2  where  fa 
stands  for  sampling  frequency,  extremely  sharp  cut-off 
slopes  may  be  obtained.  As  an  example  plot  of  cut-off  slope 
of  bilinearly  transformed  Butterworth  filter  v.s.  cut-off 
frequency  is  given  in  FI3.  34. 


C.  HEQUIHED  DATA  CARDS 


The  data  cards  required  to  use  the  program  are  given 
below . 

Card  1 ; Values  of  n,m,)c  in  312  format 
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Figure  34 


5/2 


6jT^-o{f 


- CUT-OFF  SLOPE  OF  BILINEARLY  TRANSFORMED 
BUTTERMORTH  FILTER 
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Card  2 ; Location  of  inserted  zeros,  w;  where  i = 1,2,...,in 
in  8F10.5  format 

Card  3 ; Initial  and  final  values  of  the  frequency  to  be 
used  for  the  frequency  response  plot  in  2F10.5  format 

Card  4 ; Number  of  solutions  required  in  13  format 


D.  REQUIRED  SUBROUTINES/FUNCTIONS 


In  addition  to  the  built-in  subroutines,  the  following 
IBM  source  library  subroutines  are  used. 

1.  POLRT 

2.  PLOTP 

3.  PSUB 

4.  PMPY 

5.  PADDM 


E.  DESIGN  EXAMPLE 

Suppose  we  want  to  design  a digital  filter  with  a cut 
off  slope  of  15  and  minumum  attenuation  in  the  stop  band  of 
6 0 Do. 

possible  solutions  for  various  types  of  filters  are 
given  in  table  II. 
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TABLE  II*  POSSIBLE  SOLUTIONS  OF  DESIGN  EXAMPLE 
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The  following  input  data  specifies  a STBC  filter  of 
order  7 with  in=2,  w#=1.36,  k=1,  sampling  period  of  1 sec., 
100  solution  points,  and  frequency  response  plot  from  0 to  3 
Hz. 
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Frequency  responses  of  three  of  the  possible 
of  the  design  example  are  given  in  Fig. 35. 
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Figure  35  - MAGNITUDE  RESPONSES  OF  (a)  MB  FILTER  N=8,  M=1, 
(b)  MC  FILTER  N=7,  M=1,  (c)  MB  FILTER  N=7,  M=2 
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In  previous  chapters  we  have  investigated  the  problem  of 
desiniag  digital  filters  in  the  frequency  domain,  i.e.  to 
meet  given  frequency  domain  specifications.  A filter 
designer  should  always  consider  the  transient  response 
characteristics  of  its  filter.  There  are  many  applications, 
such  as  digital  MTI  filter,  for  which  one  is  interested  in 
the  transient  responses  of  filters  that  are  specified  in  the 
frequency-domain.  Time-domain  and  frequency- domain 
characteristics  of  a filter  will  work  against  each  other. 


Filters  close  to  the  ideal  frequency  characteristic  can  be 


designed.  Filters  whose  time  characteristic  is  close  to  the 
ideal  can  also  be  designed,  but  filters  close  to  both  cannot 
[ 12]. 

In  this  chapter  time-domain  response  of  digital  filters 
will  be  discussed  and  it  will  be  shown  that  the  location  of 


the  transfer  functions  poles  has  a profound  effect  on  the 


transient  response  of  the  filter. 


B.  TRANSFER  FUNCTIONS'  POLES  AND  TRANSIENT  RESPONSE 


Given  a transfer  function  of  a digital 


factored  form 

M / \ ^ ^ (i-l.) 


filter  in 


where  z,  's  and  p,  's  are  the  zeros  and  poles,  respectively, 
of  the  filter.  In  general,  the  input  signal's  Z-transform 
is  of  the  form 




• C‘-v) 

where  g\  's  are  the  poles  of  the  input  function.  The 
response  of  the  system  may  be  obtained  from  the  transfer 
function  relationship 


( !-$) 

For  simplicity,  we  assume  that  all  of  tne  poles  of  the  Y (z) 
are  distinct,  and  making  partial  fraction  expansion  of  Y (z)  , 
we  obtain 


■/(*)=  = < 


_____ — 

2-qr  2.-pi 


_j 


Crt  -t 


The  response  of  the  system  may  be  decomposed  into  [5]  two 
parts  called  the  input  signal  mode  (ISM)  and  the  system  mode 
(SM)  as 


V,3M  U)  = 


4-q,  i-qi 


IfLfi 


C,i 

a-p. 


4.  _£ii_ 
a-pi  ^ a-p4 


Combining  eg.  (5-4)  with  eg.  (5-5)  and  eg. 


(5-6) 


Y(i)=  ko  -f  ViSM  (4)  + VsiM  (4) 


(...) 

(.-«) 


Eg.  (5-7)  indicates  that  the  response  of  any  linear  system 
to  any  input  will  contain  modes  generated  by  the  input 
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signals  poles  and  the  system  transfer  function  poles.  The 
time  domain  response  of  the  filter  may  then  be  found  by 
taking  the  inverse  Z-transform  of  eg.  (5-7) 


^ + ^ism  (.k)  -V  ^ k) 


where 

• • • -V  ler  ((\r)^ 

= e, 

Eguations  (5-9)  and  (5-10)  indicate  that  a fast  responding 
system  is  one  in  which  all  of  tae  system  transfer  function 
poles.  Pi  , are  sufficiently  smaller  than  unity  in  magnitude, 
in  order  that  the  system  mode  will  decay  to  zero  rapidly. 
On  the  other  hand,  a slowly  responding  system  is  one  in 
which  the  system  mode  decays  to  zero  very  slowly  (i. e.  at 
least  one  of  the  p is  close  to  unity  in  magnitude) . 


(r-)) 

H 


Those  poles  closest  to  the  unit  circle  will  be  called 
the  system's  dominant  poles  because  they  tend  to  dominate 
the  characteristic  of  the  resultant  transient  response. 


The  dominant  poles'  response  property  is,  as  stated  by 
Cadzow  [14]  "The  response  time  of  a linear  discrete  system 
is  directly  dependent  on  the  locations  of  the  system 
transfer  function's  dominant  poles.  Depending  on  the 
particular  response- time  reguirement  for  a given 
application,  we  then  have  to  correspondingly  locate  the 
dominant  poles  of  the  transfer  function.  A fast  responding 
system  necessitates  dominant  poles  of  magnitude  much  less 
than  one  ". 


The  poles  of  the  STBC  filter  for  the  orders  3-11  and  for 
various  values  of  m and  w^  are  given  in  Table  III. 
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Although  this  partial  fraction  expansion  nethod  helps  as 
to  understand  the  importance  of  the  poles  of  the  systems 
transfer  function  in  transient  response  analysis  of  digital 
filters,  evaluation  of  residues  of  corresponding  poles  of 
partial  fraction  expansion  is  not  a trivial  problem. 

He  believe  that  the  so  called  ' transfer  matrix  ' 
method,  which  we  are  about  to  discuss,  is  more  suitable  for 
digital  computer  simulation. 


Given  a digital  filter  weighting  sequence  h (n)  and  the 
input  sequence  x(n),  the  response  of  the  digital  filter  may 
be  obtained  by  convolution  summation 
(X. 


3(0 

. a'"' . 

Klo)  0 

ll(<)  ll(o) 


We  o 

in  matrix  form  as 

o • • • 0 - 

O • • • 0 

'iLdi- 

(5-11) 

VU-t)  • • • 

* w 

Or 

^ = fi.  . X 


G = 


is  the 

sy  ste 

ms  transfer  matrix. 

which  is  defined  as 

“ K(o) 

0 

0 • . • 

0 - 

K(<) 

W(o) 

0 . . . 

0 

(f-«) 

. W») 

lu-i)  • • 

K/o) 

To  find  the  systems  response  using  eg.  (5-13) , systems 
transfer  matrix  must  be  available.  To  obtain  the  weighing 
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sequence  (impulse  response^  of  the  system  to  form 
transfer  matrix  in  terms  of  the  transfer 
coefficients,  the  difference  equation  of  the 
solved  for  the  impulse  input,  i.e.. 


the  system 
function' s 
system  is 


K(fi)  = 0.0 


m 


Where  a;  and  b-  are  the  numerator  and  denominator 
coefficients,  respectively,  of  the  transfer  function  cf  the 
digital  filter. 

A computer  program ( FORTRAN  ) is  developed  to 
investigate  the  transient  response  of  the  digital  filters, 
using  equations  (5-13)  and  (5-15).  Program  listing  is  given 
in  Appendix  3. 

Step  responses  of  HTBC  filter  for  various  values  of  n, 
m,  and  w^  are  given  in  figures  36-39.  Figure  36  and  37 
indicate  that  increasing  the  values  of  m and  also 
increases  the  overshoot  and  settling  time,  but  doesnt  have 
any  significant  effect  on  the  rise  time.  The  rise  time  tends 
to  increase  with  increasing  order  of  the  filter. 
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C.  SUaMARY 


Transient  responses  of  digital  filters  depend  on  the 
position  of  the  poles  of  its  transfer  function.  A fast 
responding  filter  has  poles  of  magnitude  much  less  than  one. 
Modification  of  all  pole  filters  increase  settling  time, 
decrease  peak  overshpot,  and  doesnt  significantly  affect  the 
rise  time  of  the  filter’s  step  response. 
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Figure  39  - STEP  RESPONSES  OF  B,  MB,  MTBC  FILTERS ( N=5, 

a=1,  W =1*36) 
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A?PENDIX  A 


COMPUTER  PROGRAM  LISTING 


CIMENSICN  PNUMR(25  » .RNUM  I ( 2 5 ItRONUMR  <25  > ,fi  CMMI  ( 25)  » 
*C1(25)  .C<25)  ,RH25)  ,P(25) 

CnENSlCN  X(50),Z(50),I^(50)  ,Y(50) 

CI^'E^SIC^  FCCTR(5C)  ,RGOTI(5C),CaF(5C) 

CINENSICN  C^<50)  »CV(5C) ,CZ(50) 

CIMeNSlCN  RCOTDR < 50 > ,RQOTC I ( 5C ) 

CIMEASICN  CYK50) 

CIMENSICN  F>(  100)  tPV(  ICC)  fPC  (ICO)  tRCK  100) 

CIMENSICN  FN(25) ,RC(25) 

OIMENSICN  >T(25) ,YT(25),AX(25),AY(25) 

CIMENSICN  CUMMYC  10 0 ) , CUMMX ( 1 30 ) , I T8 ( 2 3 ) ,RTB < 30  ) 
CIMENSICN  CbMKAY(lCC) 

CCMPLE)  F1,R,C,C1 

f2ac  the  cfcef  cf  filter, crcer  cf  inserted  zercs  and  the 

WEIGHTING  FACTCR 

F5AC (5  ,5C1  ) N,M,K 
5C1  FORMAT<3I2) 

peac  lccaticn  cf  zeros 

FEAC(5,5CC)  (w(I),I=l,M) 

5CC  FORMAT  (8F1C. 5) 

REID  INITI/L  ANC  FINAL  VALUES  CF  THE  FRECUENCY 
no  BE  LSEC  IN  FREQUENCY  RcSFCNSE  CALCULAT ICNS ) 

PEAC (5  ,5C2  ) WBEGIN,WLAST 
5C2  FCRMAY  (2F1C.5) 

REAC  THE  NLMEEP  CF  SOLUTION  POINTS 

PEAD(5,5C:)  NPCINT 
5C«.  FCRMAT  ( 12  ) 

PEaC  sampling  period  cf  the  DIGITAL  FILTER 

READ  (5,70  T 
7C  FORMAT  (FlC. 5) 

XPCINT=NFCINT 

WCELTAMWLAST-WBEGIN  )/XPOINT 
CAcCLLATICN  CF  NUMERATOR  FCLYNCMIAL 


ICIMV=1 

Y(l)=i.C 

IF (M  .EC  .C  ) GO  TO  401 
CC  2 1 = 1 ,N 
X(1)»V«(1  )’*’»4.0 
X ( 3 ) =-2  .C^W ( I )*»2.C 


X(2)  =C.C 
X(5)=1,C 

I Sf  ) sQ 

CALL  FNFY(2,IDIMZ,Y,ICIMY,X,5) 
CC  3 J=1,ICIMZ 
Y ( J) =Z  (J  ) 

CGNTINLE 
ICIMY=  IC  IMZ 
CCNTINLE 


CCWVERSICN  CF  INDEPENDENT  VARIABLE  FROM  W TO  S 


CO  26  1 = 3 ,ICIMY,  4 
Y(I )=-l.C<Y(I) 
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26  CO^aINLE 

FIND  THE  FCCTS  OF  NUI^ER/TCR  SfU/RcO  FUNCTION 
CALL  PCLRT  I >f  ,C0 F , HNL N , RG CTR  , RCQTI  , I ER  ) 


SELECT  THE  RIGHT  HALF  PLANE  RCQTS  OF  NUMERATOR  SCUARED 
FUNCTION 


CO  20  I=1,MNUM 

IFCRCCTR (I  ).GT.0.0)  GO  TO  20 
SNUMR ( L ) =FCGTR( I ) 

RNUMKL  > = RCCTI(I  > 

L=L+1 

2C  CONTINLE 
MN=MNUN/2 
LK  = L 
LZ=L-1 

IF(LZ.EC.MN  ) GO  TO  22 
LL  = 1 .C 

CO  35  I=1,MNLM 
L)(  = L4l 
L'f=I  + l 

IF (RNUM I ( Ll ) .EQ.ROQTI ( I ) .OR .RNLM I ( LL I . EC  . ( - 1 .O^ROOTI 
» (I) ) ) GC  TC  35 
IF(LL.EC.l)  GO  TO  36 
LJ=LL-2 

IF(RNLM(LJ).EC.RCCTI(I)  .OR  . RNUM  I ( L J ) .EC.(  -1  .C=*ROOTI 
»(I) ) ) GC  TC  35 
IF(LL.EC.2  ) GO  TO  26 
LI =LL-^ 

IF  (RNUM  I ( LI ) .EQ.kCCTI  (I  ) .OF.RNLMK  LI  ) . EC  . { - 1 .0*fiGCTI 
*(!)))  GC  TC  25 
26  RNUMR( L)=RCCTR(I  ) 

RNUM I ( L ) = RCCTI ( I ) 

RNUMR(LX) =RCCTR(LY  ) 

RNLMI ( LX ) =RCCTI( LV  ) 

IF(LX  .EC.MN  ) GO  TO  2 2 
LL=LL42 
35  CONTINLE 
32  CC  22  ]=1,MN 

R(I>  =CNFL>  (RNLMR  (I  ) ,RNUM  ( I ) ) 

22  CONTINLE 

CALL  MAKFCL  (MN,R,C) 

MM  = MN-tl 

C(MN1)=CMFLX( 1.0, C.O  ) 

CC  31  I»1,MM 
RC(I >=REAL(C(I) ) 

31  CCNTINLE 
GC  TC  5 
4CJ  RC(1»=1.C 
MM=1 
MN=C 

RNUMR ( 1 ) =C.C 
RNUMI ( 1 l = C .0 

calCllaticn  cf  cenominatcr  polynomial 

5 N0NUM=2^N 

IFIM.EC.C  .ANC.K.EG.N  ) GO  TC  402 
CCNST=1.C 

IFIM.EC.C)  GC  TO  7 
CC  11  1 = 1, N 

CMLLT=(V.(I)**2.0-1.01>*=»2.0 
CONST=CMLLT+CONST 
11  CONTINLE 
/ IF(K.EC.C)  GC  TO  52 
KK=2*K 
KKK  = KK  ■»  1 
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CM(KKK  )=CCNST 
CG  12  1 = 1, KK 
C^'(I  )=C.C 
12  CGNTINLE 
GO  TC  52 

52  CM(1)=1.C 
KK  = 0 .0 
KKK=1  .C 

53  NK=N-K 
NKK1  = M<*  1 

GETAIN  Tt-E  Ct-EcYSHEV  FCL  VNCM  I A L , SCU  AR  5 IT  AND  NULTIPLY 
TEE  RESULT  EY  eUTTERk»ORTI-  SQUARED  FUNCTICN 

CALL  GEES  V (CYl ,NK) 

CALL  FEPY  (CY,NKK,CY1,NKK1,CY1,NKKU 
CALL  FEFY  < C 2 , 1 D I MC Z , C Y ,NKK  ,0^  , KKK  ) 

CETAIN  CE^C^<I^JTOR  SCUAPED  FUNCTICN 

CALL  PACCE<Z,IDIE'Z,Y,ICIMY,1.3,CZ,IDiyC2) 

CCNVERSICN  CF  INDEPENDENT  VARIABLE  FROM  W TC  S 

i DC  17  I=3,ICIMZ,A 
za)=-i.C’»z(i) 

17  CCNTINUE 
E N0NUN=IC1N 2-1 

FIND  THE  FCCTS  CF  DENOMINATOR  SQUARED  FUNCTION 

CALL  PCLRT  ( Z , CGF  , NO  NL'M  , RGGTCR  ,RCOTC  I , I ER  ) 

GC  TC  ACE 

AC2  CALL  RCCT (NCNUM, RCCTCR fPCQTCI  ) 

select  tee  right  half  PLANE  RCCTS  OF  CENCMNATGR  SQUARED 
FUNCTION 

AC3  L=1 

CO  21  I=1,NCNUM 

IF (RCCTCR ( I ) .GT  .0 .0  ) GG  TO  21 

RDNUMR ( J )=RCCTDR  (I  ) 

RDNUMI ( J ) =RCCTDI ( I ) 
j = Jf  1 

2]  CCNTINUE 
MC=NCNUN/2 
CG  22  1 = 1 , EC 

RKI  )=CNFL>(RDNUMR(  I ) ,RCNUM  (I  ) J 
23  CONTINUE 

CALL  MKFCL  (MC,R1  ,C1  ) 

NCl  = MC-»  1 

CKMCl  l=CEFLX(  l.CrC.C  I 
47  FORMAT  (2X  ,25F11. 3// ) 

CO  32  1=1, MCI 
RCKIMREAKCKI)) 

32  CONTINUE 

NCRMALI ZAT ICN  CF  NUMERATOR  POLYNOMIAL 

FACTGR=RC1(1J/RC(  1) 

CC  ICC  1=1, MM 
RC(I ) =PC (I  )*FACTCR 
IOC  CONTINUE 

WRITE (E  ,57C  ) 

CUT  PUT  SECTION 

57C  FORMAT  (•  1 ' ////) 

WRITE(E,5EC  » 

5tC  FCRMAT(2X,'  LOW-PASS  PROTOTYPE (CONT l nuGUS ) FILTER'//) 
WRITEIA  ,571) 

571  FQRMAT(2X,'  ORDER  CF  FILTER  '//) 
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V»RITE(«,eCC)  N 

ecc  F jrlM^T  (2X  , 12// > 

,572  ) 

572  FORM4T(2>,'  OROER  CF  INSERTED  ZEROS  •//) 

WRITE  (6  , ECC)  M 

WRITEIE  , 57  3 ) 

573  F0RMAT(2>,'  WEIGHTING  FACTOR  OF  TBC  FILTER  '//) 
WRITE(6,6CC)  K 

WRITECf  ,57A) 

574  FCBHAT(2>,'  LOCATICNS  QF  ZEROS  •//» 

WRITEIt  ,6C1  ) (w(  I ) ,I  =1,M 

8C1  FORMATIEFIC.S//) 

WRITE (6  , 576  ) 

576  FCRMA7(2X,'  COEFFICIENTS  GF  NUVERATOR  PCLYNCMAL 
*(  CESCENCING  ORDER  ) '//) 

WRITE(6,47)  (RC{  I ) ,I  = 1,I«N1  ) 

WRITE (6,575) 

575  FORMAT(2X,'  COEFFICIENTS  OF  CENO)'I\ATCR  POLYNOMIAL 
♦(  CESCENCING  CROER  ) '//) 

WRITE(6,47)  (RCKI  ) ,I  = 1,MD1) 

WRITE  (6  ,511  ) 

511  FORMAT  (2>  , 'ROOTS  CF  NL^-ERATGR  ' // ) 

WRI7E(6,512) 

512  FORMAT  (ex  REAL  P ART ', AX ,' IMAG  INARY  PART  •//) 

CO  52C  1 = 1, NN 

WRIT6(6,512)  RNUMR(  I ) ,RNUMI  ( I ) 

513  FORMAT (5X  ,F10.5,AX  ,F1C.5/ ) 

52C  CONTINLE 

WRIT6(6,514  ) 

514  FORMAT ( 2X  , 'ROOTS  CF  C E NCM I N ATOR ' // ) 

WRITE(6  ,5  12  ) 

CC  521  I=1,MC 

WRITE(6,513)  RDNUMR(I)  ,RDNLMI(I) 

521  CCNTINLE 

WRITE  (6  , 5 76  ) 

576  FORMATI2X,'  INITIAL  YALLE  CF  FRECUENCY  '//) 
WR1TE(6,51C)  W6EGIN 

5IC  FORMAT  (2X  ,F1C. 5// ) 

WRIT£(6,577) 

577  FCRMAT(2X,'  FINAL  VALUE  OF  FRECUENCY  •//) 

WRITE(6,51C)  WLAST 

FLuTTING  SECTION 

WRITE! 6 ,5«C) 

59C  FORMAT  ( • 1 ' J 

call  PPLCT  TG  find  FRECUENCY  RESPONSE  VALUES  CF  THE 
F ILTER 

CALL  PFLCT(RC,RC l,WBEGIN,WCELTA,NPOINT,PX,PY,MDl  ,MN1) 
CALL  FLCTF  ( FX , PY  , N FC  INT , C ) 

PRECISTCRTICN  FOR  BILINEAR  TRANSFORMATION 

SCALE=TAN  (7/2. ) 

CO  777  LK=1,MN1 
EXF=LK-1 

RN(LK)=RC(LK)/(SCALE**6XP) 

777  CCNTINLE 

CO  776  LK=1,MD1 
EXF=LK-1 

RO(LK)=flCl(LK)/(SCALE**EXF) 

776  CONTINUE 
MT=MN 
MK=M0 

CALL  ZCMN  (RN,RD,MK,MT,WBEGIN, WLAST, WCELTA,NFCINT,T) 
CALL  XFCFM(RC,RCi,l,XT,YT,MN,MC,lC0C.C,12.C) 

CALL  FFLCT(XT,YT,5C0.C,10.C,100,AX,AY,MC1,MM) 

CALL  FLCTF(AX,AY,1C0,C) 

STCP 
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C 

C 

C 

C 

C 

C 


SL6R0LTINE  CHBSV  (CY,^C) 

SCoRQLTINE  7C  FIND  THE  CCEFFICIEMS 
fclyncmial  cf  given  CPCER 


OF  CHEEYSHEV 


NC  : CRCEF  CF  CHEBYSHEV  FCLYNCMIAL 

CY  : CALCLLA7EC  COEFFICIENTS  (CESCENOING  CRCtR) 


GO  TC  2 


CIMENS  ICN  CX(50i  ,CY(5CJ ,YY(50) ,ZZ( 50  ) ,Z( 50) 
NN=NC-»  1 
CX(l)=l.C 
IF(NC.EC.C) 

CY(1)=C.C 
CY(2  > = 1.C 

IF(NC.EC.l)  GC  TC  19 
YY ( 1 ) = C .C 
YY(2)=2.C 
CO  5 I=2,NC 
II=I+1 
I 1 = 1-1 

CALL  PNFY(2,IOIMZ»YY,2,CY,I) 

ICIMZZ=ICINZ 
CC  6 J = 1 f ICIMZ 
ZZ(J)=Z(J) 

1 CONTIME 

call  F<La(Z,IDIMZ»ZZ,IClMZZ»CX,Il) 

LFGATE  FCLYNCMALS 

CO  7 J=1,I 
CX(  J ) = CY ( j I 
7 CONTINLE 

CO  8 J=  1 ,ICINZ 
CY(  J > = Z (J  I 
£ CCNTINLE 
f CONTINLE 
1C  RETURN 

2 CY(l)=CXn) 

GC  TG  H 
ENC 


SUBRCUTINE  ZCMN  ( RC  , R C 1 1 N , M , F U F2  , WOE  LT  A , NP  , ) 

C 

C SLdRCUTINE  TC  FIND  ANC  FLCT  FREQUENCY  RESPONSE  CF  DIGITAL 
C FILTERS 

C RCl  ; COEFFICIENTS  CF  OENCMNATGR  CF  CONTINUOUS  FUNCTION 
C PC  COEFFICIENTS  CF  NUMERATOR  CF  CCNTINUCLS  FUNCTION 
C N : ORDER  CF  CENOMINATCP 
C M ; CRCER  CF  NUMERATOR 
C FI  : INITIAL  FREQUENCY 
C F2  : FINAL  FREQUENCY 
C T : SAMPLING  PERIOD 
C 

CIMENSICN  ECC(25) 

DIMENSION  PC(25) ,RC1(25),ZX(25>,ZX1( 25)  ,XJ(25> ,YJ<25>  , 
»«RlOTPC  (2  5 ) ,RCOTRN  (25  ) ,RCOT  IC(25  ) .RCOTIN  ( 25  ) ,CY(25) , 
^X0AY(2f  ) , >FAYDA( 25  ) , 

♦ CX(25)  ,YFAY ( 3 ) tYPAYDA  (2  I 
CIMENSICN  ACC(iO),Ail(10»,A22(10),Bll(10),B22(10» 
CIMENSICN  XAXIS( 15C0), YAXISK  1500), YAXIS2(  15  00), 
*FM15CC  ) ,CEL(  1500  ) 

CIMENSICN  >C0F(3)  ,X2CCF(3)  , COF ( 3 ) , RGCT R ( 2 ) , RCOT ! ( 2 ) 

1 ,RC0T1  (2  ) ,RCCT2( 2 ) 

CIMENSICN  CUMKAY(IOO) 

C0U8LE  PRECISION!  W,Y 

COMMON/ FA). /W(  100  ) , Y ( ICC  ) , I CALL 

REAL  MAGH  ,MAGNH,NUMi  ,NLM2 
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CCMPLE>  C>,CYfXJ,VJ 
CCf'PLE)  MNtCENOM 
MJ=M+1 

CALL  E2>FPV(RC,M,NyM,ZX,l) 

CALL  E2XFP^(BCl»Nf^^'.‘^^ZXl,C) 

SFACT=ZX(^J)/ZX1(^J ) 

CC  <50  1 = 1, 

ZXK  I)=2X1(I)*SFACT 
9C  CCNTINLE 

V«PITE(£  ,£C) 

CCEFFICIEMS  CF  CENOyiNATOR  PCLYNOMIAL 

WRITE(£,45  KZXK  I ),I  = 1,NJ  ) 

V«RITE(£  ,£C) 

CALL  PCLR7  (ZXtCOF, Nt  FCOTRN  ,RC07IN, I ER ) 

CALL  FCLF7 (ZXi.CCF  ,N  ,RC07RC,RCC7ID, I ER ) 

CALL  CCE\<E^  (N,MX) 

CALL  FAC7CR  ( R007RN , RCG7 IN , N , N X ,A0D,A11,A22» 

CALL  FAC7CR ( R007 RC , RCC7 I C , N , MX , 322 , B 11 , ECO) 

F = F1 

CO  10  1=1, NF 
> = F’#‘7 
MAGh=l.C 
FF3  = C .C 
CO  6 0 = 1 .I'-X 
AC=ACC ( J ) 

A1  = A11  (J  ) 

A 2 =A  2 2 ( J ) 

E1=B11(J) 

64=62210) 

NLM1=(A2  + AC)*C0S(  >)+Al 
NUM2=  ( AC-A2  )<*SIN  (X  ) 

NLy  = CP  FLX  (^LM  ,NUy2  ) 

CENOMl*  (E24  1.0)*CCS( X)+B1 
CEN0y2= { 1 .-E2  )*S IN ( X ) 

DENay=CNFLX(DENQyi,CENCy2l 
MAGH=AES  (yAGF*CA3S (NLN/OENCM) ) 
FH1=(ATAN(NLN2/NUM1))<*57.2S57795  1 
PH2=(A7AN  (CENOM2/CENCM)  )*57.295  779  51 
FF 2 = FF  - ♦ ( FF1-PH2  » 

£ CGN7IME 

MAGNF  = 2C  .*ALCG10( NAGF ) 

FF  ( I ) = FF2 
FK=Ph ( I ) 

>AXI  S ( n = F 
VAXISl  ( I )=FAGH 
YAXI S2 (I ) =FAGNH 
V.(I)=XAXIS(I) 

Y(I)=YAXIS1(I) 

F=F+VDEL7A 
1C  CCN71NLE 
V»RI7E  (£  ,£C  ) 

6C  F0RMA7('l') 

Af  FGRMA7  (2X , 1QF12.3/// ) 

)»BI7E  (£  ,7C  ) 

7C  FGRMA7(2>,'  8IL1NEARLY  7RANSFGRMED  LCW-PASS 
♦(FRC7C7YFEI  CIGI7AL  FIL7ER'//) 

W3I7E(£  ,£1) 

61  FCRMA7(2X,'  C0EFFICIEN7S  OF  NUyERA70R  FCLY^CMIAL 
V«PI7E  (£  ,Af  ) (ZXd  ) , I = i ,NJ  ) 

)»RI7E  ( £ ,£2  ) 

62  FGRMA7 (2X  , 'CGEFFICIENYS  OF  OENOMINATCR  FOLYNCMIAL 
FRI7E(£  ,AE  J (ZXK  I ) ,I=1,NJ) 

V.RI7E(£  ,62  ) 

62  FGRMA7(2X,'  NUMBER  CF  CASCADED  S7AGES'//) 
V<RI7E(£,£A)  FX 
6A  F0RMA7  ( IX  , 12) 


' //  ) 
• //  ) 
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81 

83 

83 

65 

66 
6 7 

68 

9? 

13 

51 

16 

53 

1-^ 

53 


V»BITE(6  ,81) 

FORMAT  (2X, 'COEFFICIENTS  OF  CASCADED  STAGES'//) 
DC  32  1 = 1, NX 

V»RITE  ( 6 ,83  )A00(  I ) ,A11(  Ii,A22(I  ) , 81 1 ( I ) , 622  (I  ) 
FORMAT  (2>,5F12. 5//  ) 

CCNTINLE 
WRITE(6,65  ) 

FORMAT  (2X  , 'ROOTS  CF  NUMERATOR'//) 

DC  66  1 = 1 ,N 

WRITE(6,45)  POOTRNd  ),ROOTIN(I  ) 

CONTINUE 
V«RITE  (6,67) 

FQBMAT(2X,'  ROOTS  CF  C6NOMINATCR  '//) 

CG  68  1 = 1, N 

k*RITE(6,A5)  ROOTRD  ( I ) , ROOT  I D ( I ) 

CONTINUE 
WRITE (6  ,98  ) 

FORMAT(2>,'  SAMPUING  FRECUENCY  •//) 

WfiITE(6 ,99  )T 
FORMAT  (2X  ,F1C. 5//) 

N=NP 

WRITE (6  , 13  ) 

FORMAT  Cl') 

CALL  FLCTF ( XAXIS  ,YAXIS1  ,N  ,C  ) 

WRITE ( 6 ,51  ) 

FCRMAT(//^CX,'  A8S.  GAIN  V.S.  FREQUENCY') 
WRITE(6,  16  ) 

FORMAT  (•  1 ' ) 

CAUL  FUCTF(XAXIS,YAXI£2,N,C) 

WRITE  (6,53) 

FORMAT  (//ACX  ,'  GAIN  (C3  ) V.S.  FREQUENCY  ') 
WRITE  ( 6,0 
FORMAT  ( • 1 ' ) 

CALL  FLCTF(XAXIS,FF,N,C) 

WRITE(6,53  > 

FORMAT  (//ACX , ' PHASE  V.S.  FREQUENCY  ') 

RETURN 

ENC 


SU8RCUTINE  CCEVEN(N,M  ) 

SUoROUTINc  TC  FIND  THE  NUMBER  CF  CASCACEC  STAGES 
N;  CEGREE  CF  NLMERATCfi  FCLYNQMIAL 
M:  CEGREE  CF  CENOMINATCR  FCLYNCMIAL 

NN=(N/3)*2 
IF(N.Nc.NN)  GO  TC  1 
M=N/2 
GC  TC  3 
1 M=(N+1)/Z 
3 RETURN 
ENC 


SUSRCLTTNE  EZXFRM(X,M,NMM,ZX,IFACT) 

CIMENSICN  X(25),Y(25),ZX(25),XCUM(2) 

= ,YCUM(3), 21(25), Z2(25),X1(25),X2(25) 

^ ,XFACT  (35  ) ,2X1  (25  ) 

C SUBROUTINE  TC  FIND  THE  BILINEAR  Z-TRANSFCRM  OF  A GIVEN 
C FClYNCMIAL 
C X:  FCLYNCMAL 
C M:  CRCER  CF  FCLYNCMIAL 

C NMM»N-M  ( THE  DIFFERENCE  BETWEEN  THE  CRCERS  OF  NUMERATOR 
C AND  CENCMINATCR 
C ZX:  RESULTANT  FCLYNQMIAL 
C IFaCT=  1 FCB  NUMERATCR 
C IFaCT=  C FCR  CENOMINATCR 
C 

NM=M+NNM+1 
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N2=M ♦! 

CO  I 1 = 1, 

ZX(I ) = c .c 

1 CQNTIME 
>DLy(l)=-l.C 
X0UM(2  ) = 1.C 
\DLM(1)=1.C 
VDLM(2)=1.C 
IZX  = 2 

CO  2 11  = 1, N2 

CALL  PCLE>F(XDUM,2,N'N,X1,I  n 
CALL  FCLEXf(VCU>' ,2 ,1^  ,X2, 12  ) 

CALL  PNPV(Z1,IZ,X1,I1,X2,I2) 

FACT=X (II ) 

CALL  FADCf'  (Z2  ,IZ2,ZX  ,IZX,  FACT,  Z1  ,IZ) 
CC  4 K=1,IZ2 
ZX(K  ) = 22  (K  ) 

< C3NTINLE 

IZX=IZ2 

2 CCNTINLE 

IF(IFACT-EC.C)  GC  TC  5 
CALL  PCL£XF(YCU>1,2,NMM,XFACT,iNXFACT) 
CALL  f>FY(2Xl,MZ,XFACT , N X F ACT , ZX , I Z2 ) 
CC  6 1 = 1, F2 
ZX( I » = 2X1(  I ) 

6 CONTINLE 
f Ni  = M2-l 
PETURN 
ENC 


SLBRCLTINE  FACTOR ( R R ^ , P IN , V , VX , A22 , A 1 1 , AGO J 
CIMENSICN  RRN(25),RIN(25),XR(25),XI(25),YR(25),YI(25), 
» C(25»,C1(25),A00(25),A11(2£I,A22(25» 

CCNPLEX  C.Cl 

INITIALIZE  CCLNTERS 

IM  = 0 
1Y  = 1 
K = 0 
M=0 
IX=1 
NXR  = C 
NXI  = C 

IF(!^,NE.C  ) C-C  TO  18 
CC  17  I = 1 ,NX 
AOO(I)=l.C 
AIK  I ) =C  .C 
A22(I)=C.C 
17  CCNTINLE 
C-C  TC  2C 

IE  IF(M.NE.l)  GC  TO  16 
AOOd  ) = 1 ,C 
AIK  1 » =-  1 .C«RRN(  1 ) 

A22(  n=c.c 

K = 2 

NCIF=M> 

NXK=  1 
GC  TC  1C 
K CO  1 1 = 1, ICC 
1X1  = 1XK 

IFdX.GI.M  GC  "^C  7 

IF(RRN(IX).EC.RRN(IX1).AND.RIN(IX).£C. 

♦ (-l.C*FIN  ( IXI) ) I GC  TC  2 
IN=IN+1 

XR ( IN ) =-  I ,C^RRN( I X ) 

IX=I X+1 
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NXR=NXP+1 
GC  TC  1 
2 m = IV4l 

VR(IY)sRR^(IX) 

YI(IY)=FI^(IX) 

YR(IY1)=RRN(IX1) 

YI  (IYU  = RI^(  IXU 
IX=IX1-»1 
M=NI42 
^XI=^XI4l 
IY=IY1+1 
1 CQNTIME 
7 M = NXR-*NXI 
IF(NT.EC.^X) 

IF(MX.GT.N7  > GO  TO  a 
^CIF=^^-^X 
IC  = 1 

CO  4 I=1,NCIF 
IC1  = IC-H 
ACC(  I ) = 1.C 

Alim=XR(lC>+XR(ICl) 

A22(  I ) =XR  nC  )*XR  ( ICl ) 

IC=IC+1 

4 CGM1M.E 
^0IF1=^CIF41 
NCCF  = NC  IF42 
^LAST  = ^X-^  XI 

CC  5 I=NCIF1  ,,NLAST 
NCCF=NCCF4  1 
ACC( NC  IFi  ) = 1 .0 
AIKNCIF  1 ) =XR(NCGF  » 

A22(NCIF1  )=C  .0 

5 CCMIME 
N6EGIN=NLA£T+1 
NCCMF=  1 

CC  6 I=^EEGI^,MX 
NCDHPl^NCCI-F  + l 

C(l)=C^FLX(Yfi(,NCaMF),YI(NCC^'P)) 

C(2)=C^FLX(YR(NCawFl),YI(NCL^<Pl)) 

CALL  ^AKFCL(2,D,0n 

A0C( NBEG  IN  ) = 1 .0 

All{NeEGIN)  = PEAL(Cl(2)  ) 

A22(NE£G  IN  ) = REAL (Cl(  1 » ) 

NC0'^F  = NCCN  Fi  + 1 
t CGNTINLE 
GC  TC  2C 
£ NCIF=NX-NT 

IF(NXR.EC.C)  GO  TO  15 

CC  9 1 = 1, NXR 

AGQ(  I ) = 1 .C 

AIK  I ) =XR  ( n 

A22(  n=c  ,C 

K=K+  1 

« CCNTINLE 

IF(NXI.EC.C)  GO  TC  1C 
15  NXK=NX? 

CC  11  1=1, NXI 
NXK  = NXI<  + 1 
IK=I41 

Cm=CNFLX(YR(I),YI(in 
0(21  =CNPLX  (YR(IK»  ,YI  ( IK  )) 

CAlL  MAKPCL  (2,D,C1) 

ACC<  NXK  K 1 .C 

AIK  NXK  » =FEAL  (01  (2  ) ) 

A22(NXK  ) = FEAL(01(  11) 

K = K+  1 

11  CCNTINLE 
1C  CO  12  I=K,NCIF 
^XK  = NXK4  1 
ACC(NXK)=1.C 
AIKNXK  l=C  ,C 
A22(NXK  ) =C  .C 
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12  CCNTINLE 
GC  TC  2C 
3 IX=1 

CC  13  1 = 1, NXR 
AC0( I ) =1.C 
AIK  I)  = XP  n ) 

A22( I ) =C  .C 
1-  CONTINLE 

NXX=NXR41 
CC  lA  I=NXX,VX 
IX1=IX4  1 

C (K  =CRFLX  (VR(  IX  I , YI  ( IX  n 
C(2)  =CRFLX  (YRdXi  » ,YI(IX1)  ) 
CALL  MAKFCL  (2,D,C1) 

ACC(NX>  ) = 1 .0 

AIK  NX>  ) =REAL(01  (2  ) ) 

A22(NX>  I =PEAL(DK  1)  ) 
IX=IXl4l 
NXX  = NX>+  1 
14  CCNTINLE 
2C  RETURN 
END 


SLBRCLTINE  PCLEXP(XX,IOINXX,Ni,YY,IClNYY) 

SL5RGLTINE  7C  FIiNO  THE  FCV^ERS  CF  GIVEN  FCLYNQVIAL 

XX  : FCLYNCMAL 

ICI  P'XX  : CINENSICN  GF  XX 

>*i  FQViER  1C  EE  RAISED 

YY:PESULTANT  FCLYNGMIAL 

ICIPYY:  CINENSION  OF  RESULTANT  PCLYNGMIAL 

CIHENSICN  XX (25  I , YY ( 25  ) ,ZZ ( 25  ) 

IF(M.EC.i)  GC  TO  4 
ICIMYY  = 1 
YY (I > = 1 .C 

IF(M.EC.C)  GC  TO  3 
CC  1 1 = 1, N 

CALL  FNFY  ( 2 2 , 1 DI  V ZZ  , Y Y , I C IX-Y  Y , X X , I D I MX  X ) 

CO  2 J=  1 ,ICIMZZ 
YY( J l = ZZ  ( J ( 

2 CCNTINLE 

IDIMY Y=  IC  INZZ 
1 CCNTINLE 
GC  TC  2 

4 CO  5 I=1,ICINXX 
YY  ( n = XX  ( I 1 
f CCNTINLE 

ICIMYY=ICIVXX 
- RETURN 
END 


SLBRCLTINE  FPL3T ( RC  , RC 1 , W 3 EGI N , wCELT A , N PCI  NT , PX , PY , 
=NC1,NN 1 > 

OIMENSICN  RC(100»  ,RC1(1C0) 

CIMSNSICN  FX ( 1 JO) ,PY ( ICC) 

CCNFLEX  XNLM,CNUM 
h = 1*8  EG  I N 

CC  666  I=1,NPCINT 
'*  = )«  + V»CELTA 

INITIAL  IZAT ICN  CF  REAL  ANC  COMPLEX  PARTS  OF  NUYEPATGR 
ANO  OENCMNATCR 

XNUMA  = C ,C 
XNLM6=C.C 
XNLMC=C.C 


I 
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XNLMC=C.C 

C 

XCNUf<A  = C .C 
>DNUNB=C.C 
X0NUI«C»C  ,C 
XCiNUf'CsC  .C 

INITIAL IZiTICN  OF  CCLNTEPS 

^NLMC=1 
NNUM1=2 
NMiM2  = 3 
MNLM3=A 

NCNUPC»1 
MDNUM1*2 
NCNUV2*2 
NDNUN2=A 

XC=MNLNC-1 
Xl  = NNLM-l 
X2=MNLN2-1 
X3  = f'NLN3-l 


XC0=^C^LNC-l 
XCi  = NCNLM-l 
XD2*NDNLN2-1 
XCSsf^CNONS-l 

DO  6E5  J=1,1C0 

IF(MNLN3  .LE.NNl ) XNUVC=XNUVC+RC  ( MNJ'>^3  ) ^W»*X3 
IF(MNUN2.LE.MN1)  XNLNC  = XNLNC  + RC  ( MNU“2  >>fW*‘*X2 
IFIMNUM.LE.^NU  XNU  N S =XNU^' B +RC  ( MM'JM  1 M-W**X1 
IF(MNLNC.LE.MNL)  XNL N A = XNU NA  + RC ( VNUM C ) X C 

IF(MCNLP3  .LE  .MDl  ) X0NUMC=  X CNUMC >RC  1 ( M0NLM3  J*V^»^XD3 
IF(MCMI“2.LE.>'Oi)  XCNLVC  = XCNUMC  + RC1  ( VCNUM2  )=»W=*’«=X02 
IFIVCNLMl.LE.MDl)  XCNLNB=XCNLMB  + RCL  ( f'CNLM  )=»V>*=«=X01 
IF(VCNLNC.LE.MOi)  XCNLM  A = X CNUMA  + RC  1 < (« CNLVO  ) =*  k«*^X00 

INCREMENT  CCLNTERS 

MNLMC=NNLNC*4 
MNUM 1 = NNLM  1*4 
MNbM.2  = NNLM2+^ 

MNLM3  = NMN3  + A 

MCNUMO  = NCMMO+A 
MDNL'Ml=VCNLMl  + 4 
MCNUM2=MCNLM2*4 
N0NUM3*NCNLM3+A 

XC  = XC*<i  .C 
X1=X1+A.C 
X2=X2  + 4 ,C 
X3=X344.C 


XC3=XCC-»A  .C 
XC1  = XC  14A  .C 
XC2=XC2*A,C 
XC2=XC34^,C 

IF<  ( MDNLMC.GT. MOD  .ANC.  (MNLMO.GT.MNi  ) ) GO  TC  668 
66f  CCNTINLE 

CAJ.CLLATE  MAGNITUDE  CF  THE  TRANSFER  FUNCTION 

66€  XNUM1  = XNLMA-XNUMIC 
XNLM2=>NLMe-XNUMD 
XNUM  = CMFLX (XNUM1,XNUM2  ) 

CNLMDXDNLMA-XONUMC 
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C^UM2  = >C^^-^e-X0NUMD 
CNUN  = C^FL>  (C^U^'1  »C^U^2  ) 
>MAGH  = CAES (XNUM/ONLM  ) 
FX(I >=^ 

FY(I  )=2G.C>»/SLGG10(X)“^GH) 
eti  CCNTINLE 
FETUFN 
END 


SL3RCL7INE  XFORM ( X , Y , NTYPE  , XT » YT t M , N , F 1 , F2  ) 

S16RCUTINE  FCF  FREQUENCY  TRANSFORMATIONS 

NTYPE=3  ; NC  TRANSFORMATIONS  REQUIRED 
NlYFE=i  ; LCV*-FASS  TO  LCV«-PASS  TRANSFORMATION 
N1YPE  = 2 ; LOU-PASS  TO  HIGH-PASS  TRANSFORM  AT  ION 
MYPE  = 3 ; LCU-FASS  TO  BANC-PASS  TRANSFORMATION 
N1YPE=4  ; LCU-FASS  TO  EANC-STCP  TRANSFORMATION 
X : COEFFICIENTS  OF  NUMERATOR  POLYNOMIAL 
Y : COEFFICIENTS  OF  CENCMINATCR  POLYNOMIAL 
M ; ORDER  OF  NLMERATCP  POLYNOMIAL 
N : GRCEfi  OF  CENOMINATOR  POLYNOMIAL 

XI. YT  : TRANSFORMED  NUMERATOR  AND  DENOMINATOR  FCLYNGMIAlS 

CIMENSICN  X(25»,Y(25),XT(25>tYT(25l 
NM=N-M 
N l=N  + l 
Ml=M4l 

IF(NTYFE.NE.l)  GC  TO  11 

CD  1 1 = 1, Ml 

XM=Nl-I 

XT(I  » = MI  I’tCFl^^XM) 

1 CCNTINLE 
CO  2 1 = 1, M 
>N=N1-I 

YT( I » = Y (I M (F1**XN  » 

I CCNTINLE 
GO  TO  2C 

11  IF(NTYFE.NE.2»  GO  TO  2C 
XF  =M 
YF=N 

FACT  = ( F1**YF)/(F1X-»XFJ 

K=M1+1 

CO  3 1 = 1, Ml 

XM=M  l-I 

K = K — 1 

XT(II=(X(K|4(F1**xm>  impact 
' CCNTINLE 
K=N1*1 
CO  4 1 = 1, M 
XN=M-I 
K = K - 1 

YT(  I )»Y  (K  )4(F1*»XN  » 

4 CCNTINLE 
2C  RETURN 
END 


SUBROUTINE  RGCTJN.RR.RI) 

C 

C SL3R0LTINE  TO  FIND  THE  PCCTS  CF  THE  EQUATION 
C X*»N*1=0 

C N ; ORDER  CF  POLYNOMIAL 

C RF  :ARRAY  CONTAINING  REAL  PARTS  OF  CALCLLATED  PCCTS 
C RI  :ARRAY  CONTAINING  IMAGINARY  PARTS  OF  CALCULATED  ROOTS 
C 
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CINE^SIC^  FP(25»,RI(25» 

PI  -3,IA15^.2 

XN=N 

TETA=FI/>^ 

DC  1 1*1, ^ 

XI=I-1 

ARG=(  (2.C=*Xn+i)*TETA 
RR( I ) = £I^  (AFG) 


RI  (I  »*CCS  (/!RG» 


1 CCMIM.E 


RETURN 


END 


ooo  r>oo 


APPENDIX  B 


COMPUTES  PHOGHAM  LISTING 


CQMMCN  H(  ICC) ,HM(1C0 ,100) tX(lOO) ,V(100),N, P.NPULSE 

SLbRQUTINE  TC  FINC  THE  TIHE  CCMAIN  RESFGNSc 
A DIGITAL  FILTER 

N : F.5LTEF‘S  CPCER  ♦ 1 

A : NUMERATOR  COEFFICIENTS  OF  THE  DIGITAL  FILTER 
e ; CENCMINATCR  COEFFICIENTS  ( WHERE  FIRST  COEFFICIENT 
WUL  BE  NCFNALIZcD  ANC  WILL  NCT  EE  ENTERED  TC  THE 
PROGRAM 


F = l. 

N = 24 

NFULSE*2C 
WRITE(E,1C) 
1C  FORMATCl') 
CALL  INFLLS 
CALL  HMR> 
CALL  INFLT 
CALL  CCNVCL 
WRITE((,1C) 
CALL  FICT 
WSITE(<,1C) 
STOP 
ENC 


SLERCLTINE  HMTRX 

SLBRQUTINE  TC  FIND  SYSTEMS  TRANSFER  MATRIX 

COMMON  H(  100)  ,HM(  lOJ  ,100) ,X<100) ,V(100)  ,N,P,NPULSE 

NN=N+1 

NN=NFLLSE+1 

CC  1 1=1, NN 

CC  2 J = 1 ,NN 

HH(I  ,o)=C.C 

2 CONTINUE 
1 CONTINUE 

CO  3 1 = 1, NN 
CC  4 J=1,NN 
IFU.Gl.I  ) GC  TO  4 
L=I-J+1 
HM ( I ,L  ) = H ( J I 
ti  CONTINUE 

3 CONTINUE 
WRlTE(t,2C  ) 

2C  F0RMAT(2X,'  HHATRIX  '///) 

3C  FORMAT  (2X  ,2EF5. 3// ) 

CO  5 1 = 1, NN 

WRITE(E,2C)  (HM(I  ,J)  ,J  = 1,NN) 
f CONTINUE 
RETURN 
END 


SLBRCLIINE  INPUT 

SUBROUTINE  TC  FORM  THE  INFLT  PULSE  SEQUENCE 

COMMON  H ( ICO  ),HM(  100,100) ,X( IOC), V( 100 ) ,N,P,NPULSE 
NN=N+l 


1 


noo  ooo  ooo 


NM=NPLLSE*1 
CC  i I=1,^FLLSE 
>(n=p 

1 CGNTIME 
CC  2 

xn  )=0.C 

2 CCNTINLE 
VRITE(i tlC) 

1C  FCRMAT(2X,'  INPUT  VECTOR  *//» 
V«PITE(6f2C)  (X(I),I  = l,NN) 

2C  F0RMAT(2>,26F5.3//) 

RETURN 

END 


SL0RCLTINE  CCNVQL 

SiaRQUTINE  TC  PERFORM  CONVOLUTION  SUMMATION 

COMMON  l-(10J»fHM(  100,100),  X(IJO)  ,V(1 031  ,N,f 

NN=N*1 

DO  1 I =1 

V(I  )=0.C 

CC  2 J=l,^^ 

V(I)=V(I)+HN(I,J)*X(C) 

2 CCNTINLE 
1 CCNTINLE 
V*RITE(i  ,1C) 

10  F0RMAT(2X,'  RESULT  VECTOR  '//) 

V»PITE(A,2C)  (V(I  J ,1  = 1 ,NN) 

2C  FORMAT  (i>  ,26F5.3//J 
RETURN 
ENC 


SLERCLTINE  PLOT 

SLdRGUTINE  TC  PLOT  THE  OUTPUT 

COMMON  h(  100)  ,HM(  100 ,100) , X(  100)  ,V(109)  ,N,F 

CINENSICN  T (250)  ,Y(25C)  ,Z(  250 

NN=N+1 

NCELTA=1C 

NFCINT=(NFLLSE+3)*NCELTA 

CO  1 I=1,NPCINT 

XI  = I 

T(i)  =XI 

Y(I)=C.C 

Z(I )*C  .0 

1 CCNTINLE 
K = 0 

CO  2 I=1,NFCINT,NDELTA 
I FCK.GT.NN  ) GC  TO  5 
K = K+1 

CC  3 J*l,^ 

L=I+J-1 

Y(L)=V(K) 

3  CCNTINLE 

2 CCNTINLE 

5 CALL  FLCTP(T,Y,NPCINT,1) 


PLOT  INPLT  PLUSES 
K = 0 

NM=NPLL£E-1 

CO  10  I=1,NPCINT,NCELTA 
IF(K.G1  .NN  ) GC  TG  7 
K = K+1 

CC  8 J=1,A 
L=I+J-1 


,NPULSc 


,NPULSE 
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/ 


£ 

IQ 


2(LI=X  (K  I 

CQMIM.E 

CONTINLE 

CALL  FLCTF  (T,Z,NPCIN7,3» 

PETURN 

ENC 


SLBRCLTINE  IHPULS 

SLaRQLTINc  TC  FIND  IMFULSE  RESPGNSE  OF  THE  FILTER. 

( iCLVING  CIFFERENCE  ECLATIONS  I 

C3MMCN  H(ICC) .HM(ICO.IOO) f X(100»»V(lC0)tN,P,NPULSE 
CIMEhSlCN  A(l6o»,e(lic> 

CI^<E^S1C^  CLMMY(IOO) 

CIMENSICN  C(IOO) 

REAO (5 ,2C  ) H 

NN*N+1 

NNN=NN-»2C 


IMTIALIZATICN  CF  VECTCFS 


CC  35  1 = 1, ^^N 
C(I)=C. 

H(I>»C. 

A(I)=C,  ^ 

en  »=c. 

35  CONTINLE 
C( 20) =1. 

SUFT  ORIGIN  TC  20 
N0RDER  = f‘42C 

fiEACI5,2C)(A(I », I=20,NCRCER ) 

REA0(  5 ,2C  ME(I  ) ,I  =20  .NCRCEF  ) 

FORMAT  ( 12  ) 

FORMAT  (8F1C. 5) 

V»RITE(^,1)  M 

FGRMAT(2>,'  ORDER  OF  THE  FILTER  IS  :',I2///> 
WRITE(G  ,2) 

FORMAT  (2>  , 'NLMERATCR  COEFFICIENTS  '///) 

WRITE  (6, 2 )(A(n,i=20,NCRCER) 


2C 

2C 


6C 

AC 


1C 


WRITE(A  ,A) 

FORMAT(2>,  'CENOMINATCR  COEFFICIENTS 
WRITE(<,2  ) (Ed)  , I=20,NCRCEP  ) 

CO  AC  II=2C,NNN 
FACT1  = C . 

NK=I 1*1 

CO  5C  IL=2C,NCRDER 
NK=NK- 1 

FACT1=FACT14A(IL)*C(NK) 

CONTINLE 

FACT2=C. 

NK  = II 

CO  6C  IL*21,NORDER 
NK=NK-1 

FACT2=FACT24e(IL)*H(NKJ 

CONTINLE 

H(II » =FACT1-FACT2 

CONTINLE 

WRITE  (E  ,2  1 ) 

FORMAT  (2X  , 'WEIGHTING  SECOENCE'///) 
WRITE  ( E , 1C  MH(  I ) , I = 2C,NNN  > 
FCRMAT(2X,12F10.5//) 
F0RMAT(2X,12F1C.3) 

K K K = 1 

CO  55  KK*2C,NCRDEP 
A(KKK  ) =A  <KK ) 

E(KKK  )>E  (KK  ) 

KKKsKKW^  1 


'///  ) 


100 


55  CONTINLE 
K K 3 1 

CO  65  I=2C,NNN 
CU^<MY  ( KK  )»M  1 1 
KK=KK+1 
65  CCNTINLE 
CO  75  1=1, 
h(IJ=CLNf-Y(I) 

75  CONTINLE 
>CUM  = e (2  ) 

B(3) =e  m 
e(i)=xcL>< 

PEAC (5  ,1CC  > F1,F2 
, R5A0(5,1CC)  T 

IOC  F0«MAT(2F1C.5) 

NF=1CC 

>NP=1CC. 

FCELTA=  (F2-F1  l/XNP 

CALL^ZCf-NlA  ,e  ,M,M,F1  ,F2,FDELTA,NP,T) 
ENC 
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